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We present a framework, based on the null-surface formulation of general relativity, for discussing 
the dynamics of Fermat potentials for gravitational lensing in a generic situation without approxi- 
mations of any kind. Additionally, we derive two lens equations: one for the case of thick compact 
lenses and the other one for lensing by gravitational waves. These equations in principle generalize 
the astrophysical scheme for lensing by removing the thin-lens approximation while retaining the 
' weak fields. 

o 
o 

(N I I. INTRODUCTION 

In the astrophysical approach to lensing |l|,||, background light sources are considered to lie on a plane at a 
distance D s from the observer, and the deflector or lens is described by a mass distribution on a plane at distance 
Di from the observer. Lightrays from the source in the source plane travel along straight null geodesic paths on a 
flat or cosmological background until they reach the lens plane, at which point they suffer a sharp bend in direction, 
subsequently traveling again on a straight line towards the observer. The direction at which a lightray reaches the 
observer determines the angular location of the image on the observer's celestial sphere. The celestial sphere can 
\Q 1 be identified with the lens plane (which is sometimes referred to as the image plane). A fundamental element of 
t— I , gravitational lensing is a lens equation of the form: 

o ■ 

o ti = ^e-A.a(o a) 

(N : 1 

where r) is a position vector on the source plane representing the source's location, £ is a position vector on the lens 
plane representing the location where the lightray emitted by the source pierces the lens plane, Di s is the distance 
O" 1, between the lens plane and the source plane and a is the bending angle that the lightray suffers at the lens plane, 
5^ I which is determined by the distribution of mass on the lens plane and is thus a function on the lens plane. The 
bending angle a can be obtained from a deflection potential ip(£) via 
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and the deflection potential tp in turn is obtained from the mass distribution £(£) by means of a 2-dimensional Poisson 
equation 

A^(0 = (3) 

where is the 2-dimensional Laplacian on the lens plane. 

A very elegant scheme that leads to the lens equation is based on Fermat 's principle, according to which light takes 
the path of least travel time between any two fixed points. In order to use this principle, a set of paths between the 
fixed source at rj and the observer are considered, all of which consist of two straight segments, but which are allowed 
to pierce the lens plane at arbitrary points The value of £ thus becomes a label for the path. The travel time along 
such paths is thus a function of £, which must be minimized with respect to The travel time along such paths is a 
Fermat potential 0(£,?7), in general depending on both the source point as well as the point on the lens plane. The 
Fermat potential has the form 
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D s D t 



(4) 



where the first term represents the geometric length of the broken path with respect to the length of a straight path 
between the source and the observer, and the second term represents the gravitational time delay suffered by the 
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lightray as a consequence of the presence of mass at the lens plane. The minimization of the travel time with respect 
to the path requires 

= 0. (5) 

This equation provides the lens map £ — > i~i from the lens plane to the source plane, i.e., Eq (Q). The Fermat 
potential is affected by the presence of a mass distribution £(£) on the lens plane essentially because of Eq. (^). As 
a consequence of Eq. (|J), one can see that the Fermat potential also satisfies a Poisson-type equation 

A € <MI,r,) = -S(0, (6) 

where S(£) differs from £(£) by an additive constant. This is the dynamical equation for the Fermat potential in the 
astrophysical approach to lensing. It encodes the information of the curvature of the spacetime due to the presence 
of mass on the lens plane, and it thus represents the role of the Einstein equations in lensing. One can see that the 
Fermat-potential formulation of lensing is fully tailored to the approximations in use in the astrophysical approach to 
lensing, namely: thin lenses, weak fields and small angles. 

We are interested in a generalization of the Fermat-potential approach to lensing for a generic situation with no 
approximations. The first step to this goal is to construct a Fermat potential for a generic situation, such that a lens 
mapping from the observer's celestial sphere to the lightsource locations results from a variational principle applied 
to it. We refer to this phase of the project as the kinematics of the Fermat potential. The second step is to find a 
"field equation" for the Fermat potential, namely: an equation that relates the Fermat potential to the sources of 
gravitational field. We regard this as the dynamics of the Fermat potential. 

In a previous work j| we have described a setting for an approach to the kinematics of Fermat potentials using 
both a generalization of a Fermat potential and a generalized Fermat principle consistent with a generic situation, 
which we briefly summarize in the following. 

The aim of the scheme is to construct a function on spacetime that, by extremization, leads to an expression of 
the observer's past lightcone. The observer's past lightcone can be parametrized in terms of the celestial sphere (the 
null directions at the apex of the lightcone) and an affine length along the null rays that rule the lightcone. Such a 
parametrization leads directly to a lens mapping from the celestial sphere out to lightsources, and a time-of-arrival 
equation for the light signal reaching the observer. Thus, our immediate aim is to provide a scheme by which the 
observer's past lightcone can be constructed. Since lightcones are level surfaces of solutions to the eikonal equation, 
such solutions are our starting point. 

Consider a lensing spacetime with coordinates x a and metric g a b- Suppose we have a (sufficiently generic) 2- 
parameter family of solutions Z(x a ,CiC) of the eikonal equation g ab Z, a Z,b— (with ,„= d/dx a ), such that the 
surfaces of constant value of Z(x a , £, Q foliate the spacetime for every fixed value of (£, C)- Define 

G(x%,x a ,(, C) = Z(x a , C, C) " Z(x a X, (7) 

Then, for every fixed value of (C, CX t ne surface defined by the points x a that satisfy 

GK,z°,C,C) = (8) 

is null and contains the point Xq, at which the vector £ a = g ab G,b sweeps out the sphere of null directions as (C, C) 
is allowed to vary. The parameters (£, £) are thus coordinates on the observer's celestial sphere, and can be taken as 
complex stereographic coordinates. Solving (^J) for x° = t we have 

t = T(x a ,x\(X), (9) 

namely, the value of the coordinate time of the point of intersection of the family of null surfaces G = with the 
worldline (x l , t) of a static source at the spatial position x l , £ = 1,2, 3. (This can be reformulated to include a moving 
source.) Varying ((, £), different null surfaces in the family G = are considered, all of which contain the observer Xq 
and intersect the worldline of a potential source at x l at time T(xq, x 1 , £, C)- This time cannot be interpreted in terms 
of curves joining the points Xq and (x\ T) because there is an infinity of curves lying on the surface G — joining the 
two points for which the coordinate time of the end point is t: there is no way to pick a preferred curve on the null 
surface. In fact, among all possible curves joining the two points on a given null surface G = 0, there may even not 
be one that is geodesic, and there are certainly many that become spacelike at least one point. Thus, even though 
T(xq, x % , C, C) does represent a time function in a sense, we have no reason to identify it with the time of arrival of a 
light signal. However, T(xg, x l , (, () is the time at which a wavefront of the (C,C)~ nuu surface passes by the source 
at x l , so we can think that instead of looking at travel times along null paths we are in effect looking at travel times 
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of wavefronts from the observer to the source. Both points of view are equivalent in the end, because ultimately the 
values of (C, C) f° r which the intersection time T(xq, x\ £, is stationary are null surfaces for which there exists a 
null geodesic that joins with (x l ,T) ||,||. 

Our recipe is: to extremize the function T(xq, x l , £, £) with respect to variations of (C, C) keeping x l fixed, and 
interpret the extremum T e = T(xq, x l , C e , Ce) &s the time of arrival of a light signal from x l to x l , giving an image 
in the direction (( e ,( e ). These choices of (C, C) \e&d to null geodesies connecting the observer with the source. It 
follows that the function T(xq, x l , () can be viewed as a generalized Fermat potential. On the other hand, since 
T(xq, x l , £, C), in general will be given implicitly via (||), then extremizing T(xq, x 1 , C, C) is equivalent to imposing the 
conditions 

G, c | a . = G, f | a . = (10) 

which arc, in fact, the conditions for the envelope of the level surfaces of G. This results from implicit differentiation 
d/dC\ x i of (@) where t is assumed to be given by (||). One obtains 

G, c I*- + G, t T )f \ xi =0. (11) 

Since G,t must be non-zero (otherwise it could not define T implicitly), then T,^ \ x * = G,f \ x a = 0. Likewise with 
T,£ \ x i . We refer to G(xq, x a , (, £) as an "implicit" generalized Fermat potential, making a slight abuse of terminology. 

In fact, G — plus the envelope condition for G as a joint set of equations produces for us two things. 1.) A 
lens equation: an assignment of a value of (C, C) on the celestial sphere for each source location x l , through light 
signals. 2.) A time of arrival of the light signal. Thus, the envelope conditions ( |l0| ) on a level surface of G embody 
our generalized Fermat's principle: a version in terms of surfaces, rather than paths. 

Up to this point, we have assumed that the Fermat potential G is somehow known, and satisfies the cikonal equation. 
Thus we have only considered the kinematics of such implicit Fermat potentials: we have not considered how a Fermat 
potential is affected by the presence of matter in the spacetime (the dynamics) . The matter (the deflector) is, of course, 
hidden away in the metric g a b that determines the eikonal equation, so the eikonal constitutes part of the dynamics, 
but the field equations for the metric itself are the Einstein equations R a b — \gabR — Tab- In principle, one could first 
solve for the metric, and then use the metric in the eikonal to obtain the Fermat potential G. There is no general 
solution to the Einstein equations in terms of data or matter sources, because of their nonlinearity, therefore this 
route does not have the potential to provide us with a one-step dynamical equation for the Fermat potential driven 
directly by the matter source, which would be analogous to Eq. ([|). 

Notwithstanding, there is an approach to general relativity which has the potential to directly lead to dynamical 
equations for implicit Fermat potentials |5|-|§]. We describe the resulting dynamics of implicit Fermat potentials 
through such an approach in Section 0. An application of the dynamics to the case of weak fields is worked out 



in Section III. The main results of the application to weak fields are two lens equations: one for the case of thick 
compact objects and the other one for the case of lensing by gravitational waves. Whether these equations can be of 
use in astrophysical situations remains to be seen, but our point is that they can be written down in principle. 



II. DYNAMICS OF THE IMPLICIT FERMAT POTENTIAL 

The formulation of general relativity via null surfaces |^,^] provides a setting for discussing the dynamics of implicit 
Fermat potentials in the sense of the previous Section. This formulation is described in complete detail in |||| . Here 
we provide a (necessarily brief) summary of the main aspects of the formulation, hoping that the unfamiliar reader is 
able to complement this Section with ||||]. 

In the null-surface approach, the aim is to rewrite the Einstein equations in terms of variables that are to be thought 
of as more fundamental than the metric g a b, which can then be generated from them. The fundamental objects of 
the formulation are two functions of six variables, Z(x a , £, () and f2(x°, £, Q. The variables x a represent points in a 
four-dimensional spacetime with a Lorentzian metric to be determined. The variables (£, £) represent a sphere fiber at 
every point x a ; eventually (C, C) acquires the meaning of the sphere of null directions at each spacetime point. The first 
function, Z, whose level surfaces are null for all values of (£, £) encodes the conformal structure of the spacetime. That 
is to say that if the metric of the spacetime is known, then Z(x a , (, £) satisfies the eikonal equation g ab Z, a Z,b = for 
all values of (C, C)- The second function, f2, represents the scale (or conformal) factor of the spacetime with conformal 
structure given by Z (all conformally related spacetimes have the same conformal structure; the scale factor breaks 
the conformal invariance). From the two variables, Z and f2, a Lorentzian metric can easily be constructed. Field 
equations, equivalent to the Einstein equations can be imposed on them. 
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These field equations for Z and are in general quite complicated. However, for the weak field applications that 
we will treat they become quite simple and are easily reduced to quadratures. We will first present the full equations 
before giving the linearization. 

It is perhaps easiest if the equations are expressed in two steps. We first introduce several auxiliary variables 
obtained from derivatives of Z(x a ,(,C)> wc define four variables 9 1 (x a ,C,, £) which can be considered as coordinates 
given intrinsically by the families of null surfaces by 

0° = Z = u, (12a) 

9+ = dZ = w, (12b) 

9~ = dZ = 57, (12c) 

9 1 = mz = R. (12d) 

Using these relations we can go back and forth between x a and 9 % . Further, it is assumed that these relations can be 
inverted locally so that we have 

x a = x a (0\t,t). (13) 
In addition to 9 l , we define another fundamental variable (and its complex conjugate) by 

A(0\C,C) = 9 2 Z, (14) 



where the x a on the right-hand side have been replaced by the 9 l via Eq. ([L3j). 

For the second step, the Einstein equations are written in terms of A(9 l , £, and Q{6 % , C, C) as 

n,u = Q(A)n + TQ 3 , (15a) 

m = -w(A)n, (i5b) 

M (A) = 0, (15c) 
where we have used the notation ,,; = ^ and in particular 

The quantity T encodes the information in the stress-energy tensor T a ;,, while (Q,W, M) are explicit functions of A 
and its derivatives: 

Q = — [A, ii A,n +^z(q,i f ~ 9,n ] , (17a) 



q 



M = 3(A,i ) - 2A,_ -{W + 8 In q\k,x , (17b) 
W=fA, + +ii(A )1 ) + iA a X > _+iA, 1 g(S a )-igin ? -iA > ifln« > ) (l-\K, 1 H,x\ , (17c) 



with 

g = l-A,iS,i. (18) 
The operator 3 is the usual 9 with x a held constant; but operating on functions f(9 1 ,(i it takes the form 

8/(0*, C, C) = S/Ie« + 8 + U +A/,+ +( 1 A,1 A J+ T 7 ) /,i (19) 

v 1 — A,i A,i / 

with 

J = -29+ + 0-A, o +5A| ei + 9 1 A, + +AA,_ (20) 
For the purposes of applying the 5 operator, A and f2 have spin weights 2 and 0. 
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Eq. (15a) is a direct translation of the Einstein equation, G a b = T a {>, while Eqs. ( J15q ) and ( |l5q ), referred to as the 
metricity conditions, are the requirements that a unique metric can be constructed from Z and CI. 

The field equations couple A and f2 in a complicated non-linear manner that makes finding exact solutions an 
almost impossible task. Nevertheless, as we will shortly see, the linearization becomes quite simple. In principle, the 
procedure for finding the Fermat potential would be to solve Eqs. ( |l5| ) for A = A(9\CX) = A(Z, 9Z, (3Z, S9Z, £, Q 
and then return to Eq. (p^), written in the form 

d 2 z = A(z,BZ,dz,mz,(,0 (21) 

so that it becomes a dynamical equation for Z(x a , C, C) with a known right-hand side. Actually, since A is a complex 
function, (El]) is a pair of over-determined second-order PDEs for Z in the two variables (£, C). From general consid- 
erations [p[7the solution space is four-dimensional; the solutions have four constants of integration, x a , that become 
the space time coordinates. Eq. (^lj) is our dynamical equation for the Fermat potential. In this regard, we point out 



that Eq. (15c) is precisely the vanishing of the generalized Wunschmann invariant that is central to the problem of 
the classification of this type of overdetermined systems of PDE's |^.|uj- 

Some effort has been devoted to the issue of rewriting the field equations in terms of Z and f2 by eliminating A via 
( pl| ) and (|l2|), with partial success To do this exactly appears to be a formidable task. Nevertheless, in the case 
of weak fields, the linearization of the field equations allows us to combine the content of the field equations (|15|) and 



of fl2l| ) into a single dynamical equation for Z in terms of given data We describe, in the following section, this 



application of our program to weak fields. 



III. LENS EQUATIONS IN THE CASE OF WEAK FIELDS 

The field equations ( |l5| ) have input of two types. In the first place, there is the matter and stresses, encoded into 
the symbol T. Additionally, there is also gravitational radiation at the boundaries of the spacetime that must be 
considered free input. This is because the field equations are equivalent to the characteristic problem for the Einstein 
equations, in which case the gravitational radiation incoming into the spacetime must be specified PJTl[|, usually 
in the form of the complex shear er(it, £, £)of a null congruence of the Bondi type evaluated at the boundary of the 
spacetime, J~ . This is a complex function of three arguments, the local Bondi coordinates on J~ . 

We concentrate on the case in which all input sources or data are small, namely, there exists a smallness parameter 
e by which the magnitude of the source and data can be measured, and such that its powers can be neglected. This 
will lead to a linearized version of the dynamics of the Fermat potentials discussed in the previous section. However, 
within this regime there are two distinct main cases of interest. One is the case where the spacetime is empty of matter 
and the gravitational radiation is small, that is: T — and |er| ~ e. We refer to this case as the gravitational-radiation 
case. The other case is such that the matter source is weak but finite, and the gravitational radiation is negligible, 
namely: \T\ ~ e and \o~\ ~ e k with k > 2. This case will be referred to as the near- Newtonian case. 



A. The gravitational-radiation case 



We assume that we have a purely radiative spacetime, in which the stress-energy tensor vanishes and the gravita- 
tional radiation is small (T — and \a\ ~ e). The spacetime is then flat up to a small perturbation. In flat spacetime, 
null planes in all possible directions constitute a 2-parameter family of solutions to the eikonal equation, leading to 
Z flat (x a ,C,C) =x a L where 

e a = 1 (i + CC, C + C, -i(C - 0, CC - i) (22) 

V2(l + CO 

is a future directed null vector pointing in the (£, Q-direction. We may expect to have then 

Z(x a ,(,() = x a e a + Z(x a ,(X) (23) 

with \Z\ ~ e. Since d 2 £ a = 0, then the variable A in the field equations ( |l5| ) is of first order in e, and products of A 
and its derivatives can be neglected. Therefore the field equations ( [l5|) can be linearized with respect to A around 
A = 0. They become 
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0,11 - 0, (24a) 
3ft = ~ f A !+ +~8(A,i )) Q, (24b) 



2 V 2 

2A,_ -S(A,i ) = 0. (24c) 



By Eq. (24a) with boundary conditions cons istent with an asymptotically flat space-time (in which Q — > 1 along null 



directions), we must have fi = 1. Thus (24b) becomes 

2A,+ +3(A,i ) = (25) 

Equations ( p4oj ) and ( |25] ) can be manipulated through algebra, differentiation and integration (with boundary condi- 
tions consistent with an asymptotically flat spacetime), to yield || 

if A = d 2 a(8°, C, C) + 3 2 <x(f?°, C, 0- (26) 
Using ( |l2| ) and (|2l|), this equation can be rewritten as an equation for Z: 

t<5 2 Z = Wa(x a £ a ,(, + &v(x a l a , (, 0- (27) 



Notice that in (|27j ) the operator 3 acts at fixed x a , so that all three arguments in a(x a £ a , £, £) contribute to the eth- 
derivative. In order to obtain this equation from (^6|), <7(0°,£, £) is approximated by a{x a £ a , (,,(,), an approximation 
that is justified because a is a first-order quantity, therefore first-order terms in its arguments can be neglected. Since 
Z leads to the Fermat potential G in a trivial manner, we consider this the dynamical equation for the implicit Fermat 

potential in the presence of weak gravitational radiation. The operator S S 2 is the double Laplacian on the sphere 
(£, C)i whereas the right-hand side acts as a known source. This is the analog of Eq. (|^) for lensing by gravitational 
waves without the approximations of thin lenses or small angles. 

There exists a simple Green function [|l2| for the double Laplacian <3 <3 2 , given by 

F(cxx'X') = ^e'j a H£' a n, (28) 

which can be used to express the general solution to ( |27| ) in closed form in terms of the source: 

z(x a , c, c) = x a £ a + ±- f (H' 2 a(*x, C, C') + 3'M* a C C, C))e'J a HCHds 2 '. (29) 

This expression for Z can now be used to construct the implicit Fermat potential G = Z(x a , £, C) — Z(x$, £, (). The 
two envelope conditions 9G = 5G = then take the form 

o = (x a - x a )m a 

+g' 2 [a(.* a C C, CO - o{x%e a , C, c"')]) WJ a Hi'J a )}ds 2 '. (30) 

By our discussion in Section D, Equations (30) and Z(x a XX) ~ Z{xq,C,X) — combined are equivalent to the lens 
mapping and time delay for the case of lensing by weak gravitational waves without the approximations of thin lenses 
or small angles. Since there is no lens plane in this case, the lens equation is interpreted as a mapping from the 
observer's celestial sphere to the source's spatial location. The lens mapping is understood in the following sense. We 
assume a spin-2 function a(u, £, £) is given, perhaps constructed with some model in mind. This function represents 
the gravitational radiation that reaches the boundary of an asymptotically flat vacuum spacetime. We can specify 
our coordinates by setting the origin on the observer's worldline, so Xq — (r, 0,0,0). The point x a — (t,x,y,z) that 
satisfies both equations Z(x a , (, () — Z(xq, CX) — and ( |30| ) then lies on the observer's lightcone, and is reached by a 
null geodesic in the direction (£, £) at the observer. These are three equations for the four variables x a , so one of them 
will be free. As usual, we assume that this free coordinate represents a distance to the source, or can be obtained 
from a measured distance to the source. The remaining two spatial coordinates then specify the angular location of 
the source, and the time coordinate represents the travel time of the light signal from the source to the observer. 
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It may not be a simple matter to solve for the two angular coordinates and time from Z(x a , £, Q — Z(xq, C, C) = 
and (p0[), in order to express the lens mapping in closed form even if a function a is given explicitly. In fact, it is 
not even clear that the lens mapping can be expressed in closed form. But the point that we wish to make is that 
in principle, the lens mapping is there. Work on an explicit example of this lens mapping is in progress and will be 
reported elsewhere. 

What would be the use of lensing by gravitational radiation? Lensing by extremely compact objects provides a 
means to study the mass structure of the deflector. Likewise, we can think of using lensing by gravitational waves 
to infer information about the gravitational waves themselves. In principle, lensing events are natural detectors of 
gravitational waves. In practice, however, these lensing events are expected to lie below the observational limit, and 
be overshadowed by the magnitude of lensing effects by compact objects. 



B. The near-Newtonian case 



We turn our attention now to the opposite case, namely: there are sources of weak strength but there is little or no 
incoming gravitational radiation. This would include, for instance, the case of static spacetimes with weak (isolated) 
matter sources, or any matter sources that are static in first approximation. We thus have, by assumption, \T\ ~ e 
and \a\ ~ e k with k > 2. The source T is a function of six variables obtained from the stress energy tensor T ah by 



T(0 l , £, C) = T ab (x a )£ a £ b , 



(31) 



where 



x a = (29 a + 9 1 )£ a + . 



(32) 



where £ a is given by Q22J). In this case, the Einstein equation ( 15a ) splits into two hierarchical equations, one for the 
zeroth-order term Oq and the other one for the first-order term O: 



(fio),u = 0, 
S2,n = Tfta. 



(33a) 
(33b) 



It has been shown |jj that in the asymptotically flat regime, one can take Oo = 1 consistently with Eq. ( pL5b| ) up to 
first order. Thus we take Oq = 1 and the field equations (|l5|) become 



0,n 



T. 



30 = - I A. 



-ii(A,i 



2A,_-8(A,i) = 0. 



(34a) 
(34b) 
(34c) 



where we have used the fact that also in this case the spacetime can be represented as a small perturbation off flat 
space, in which case A is of first-order in e. By virtue of Eq. (34a), the conformal factor O can be obtained from the 
source by quadratures: 



= 1- 



OO POO 



Tde^'de 1 



(35) 



where the upper limit of the integral has been chosen at 9 1 = o o, the null b ound ary o f th e asymptotically flat 
spacetime. Thus O can be considered as a given source of Eq. ( |34b| ) . Equations (34b) and (34c) can be manipulated 
through algebra, differentiation and integration |(|, to yield 



A 



0" 



39 



(90), _ d6 u + 39 2 / (90),+ d6 u - 299(20 - 990) dff 



qO/ 



(36) 



The upper limit 0° = oo represents the point i° at the null boundary of an asymptotically flat spacetime (timelike 
infinity) . 

Using (|l2[) and (J2|), Eq. (Bq) can be rewritten as an equation for Z: 



9 2 9 2 Z 



x a £ a 



39 2 



(90), _ dO 1 ' + 39 2 



x a m£ a _ 



(90),+ d6 u - 299(20 - 990) dff 



Or 



(37) 
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This is a Poisson-type equation, since the operator in the left-hand side is the double Laplacian on the sphere (£, £), 
and the right-hand side is a given source, by virtue of (j35|). Again, since Z leads to the Fermat potential G in a trivial 
manner, we consider this the dynamical equation for the Fermat potential. This equation thus plays a role analogous 
to (||) and generalizes Eq. ([|) by removing the approximations of thin lenses and small angles. 

We can now proceed to write down the lens mapping explicitly. We start by writing down the solution of (|3^ ) in 
terms of the source in the right-hand side, by means of the Green function F((, (, (' , (') given by Eq. (^8|): 

Z(x a , C, C) - x a l a + S(x a X\ (Yj a HWdS 2 ', (38) 

with 

pOG / p OO pOO \ 

S(x a ,C,0= (39 / _ (3ft),_ dO 1 ' + 33 2 / _ (8fi),+ dd 1 ' - 253(2£1 - SSfi) ) dO Ql . (39) 

We can now construct the Fermat potential G = Z(x a , C, C) ~ Z(xq, C, C) an d write down the two envelope conditions 
3G = SG = 0, which then take the form 

= (x a - x a )m a + ^J s2 S(x a ,C, CWJ a WJ a )\dS 2 '. (40) 

This equation, with Z(x a , £, C) — Z(xq, (, () — 0, is a generalization of the astrophysical lens equation and time delay 
that removes the thin-lens and small-angle approximations. As in the previous subsection, there is no lens plane, and 
the mapping takes points on the observer's celestial sphere into the source's spatial locations. To be more specific, 
we can choose the coordinates so that the observer lies at the origin = (t, 0, 0, 0). The stress-energy tensor T a b of 
a given matter configuration on a flat background is assumed to be given, perhaps with some model in mind. The 
quantity T — T a b£ a £ b is constructed, and is used as the source of the quadrature for the quantity Q via Eq. (|35|). The 
quantity f2 thus obtained is used to obtain S according to (|39|). The combined equations Z(x a , £, £) — Z(xq, £, Q = 
and (f5J provide three equations for the four variables x a . The values of x a that satisfy the three equations lie on the 
observer's lightcone. One of the four variables x a remains free, and is assumed to be given in terms of some observed 
physical distance. The remaining three variables represent the angular location of the source and the time of arrival 
of the light signal. 

It may be complicated to write down the lens mapping in closed form, even if T a {, is given explicitly. However, one 
can envision a physical situation that is close to a thin lens, in which case we would expect this scheme to return a 
kind of "post-thin" lens equation. In other words, it is worth considering the possibility of obtaining corrections to 
the thin-lens approach by means of this scheme. The effect of the thickness of the lens has only rarely been considered 
m the literature Jlf-[l7|]. 



IV. REMARKS AND OUTLOOK 



By setting up a framework within the null-surface formulation to discuss the dynamics of Fermat potentials for 
gravitational lensing in full generality, we have achieved two main goals, as argued in the following. 

In the first place, we have provided two lens equations that generalize the astrophysical approach to lensing in two 
ways: 1.) both lens equations remove the approximations of thin lenses and small angles; and 2.) the lens equations 
include the case of lensing by gravitational waves. Even though the exceeding complexity of the case of non- vanishing 
matter sources - Eq. ( f!o| ) - puts its applicability into question, the point of our work is to show that, in principle, the 
thin- lens approximation can be removed in a meaningful manner. On the other hand, the obvious simplicity of the 
case of vanishing matter sources - Eq. ( |30|) - opens a door to the phenomenon of lensing by gravitational waves, which 
has occasionally been treated before, within the framework of the thin-lens approximation Jl8| and in a statistical 
framework []l9| . In principle, both cases could be combined in a straightforward manner, simply by adding their 
respective Z functions, therefore the spectrum of many possible applications to weak fields is opened in this work. 

One might question the necessity of the use of any type of Fermat 's principle when one moves away from the 
thin-lens approximation on the basis that the direct integration of the null geodesies on a linearized perturbation off 
flat space would clearly yield a lens map. In |L6]], for instance, the equations for the null geodesies of a linearized 
perturbation off flat space are written down explicitly and integrated to obtain the bending angle of a moving thin 
lens, whereas in jlTj the direct integration of the null geodesies of a perturbation off an isotropic cosmology is used 
to derive the cosmological thin-lens equation without the use of broken paths. Still, so far as we are aware, very 
rarely if ever has the direct integration of the null geodesies led to a scheme in any way comparable with the thin-lens 



8 



scheme in practicality. Ideally, one would like to have a correction to the thin-lens scheme that represents the effects 
of the width of the lens and an estimate of the regime in which such corrections cannot be neglected. One can easily 
argue that the approach that we present here may have the potential to make a contribution towards such a goal 
and deserves to be pursued. The usefulness of our approach must first be tested in particular cases where results are 
already available. In this respect, the thin-lens scheme has been easily reproduced in our approach p], including the 
case of moving lens planes. It would be of great interest to be able to relate our lens equations to existing results 
involving lensing by gravitational waves, such as jl9[l among others. The reader should be aware that, because of the 
statistical nature of the problem considered in |19f| , our equation ([50]) cannot be compared directly with results in 
that article. 

One should keep in mind, as well, that the treatment in this work is limited in at least two senses. We have 
only considered weak fields as perturbations off flat space. Therefore, the application to isotropic cosmologies is not 
covered here. There are reasons to think that the scheme developed here can be adapted to isotropic cosmologies 
with little effort. Work on this line is in progress and will be reported elsewhere. Secondly, we have only considered 
a coordinate-based treatment of lensing, as opposed to an optical-distance treatment. By virtue of this we have 
stayed away from a great technical and practical difficulty in observational lensing, which is the use of a meaningful 
observable distance to the source. We have bypassed this difficulty by simply using abstract coordinates. The use of 
an observable coordinate distance to the source would greatly complicate the lens equations obtained here in form, 
but it would not affect their underlying content. 

In second place, we have here demonstrated that the null-surface formulation of general relativity has the potential 
for physical applications of interest, in spite of its apparent level of abstraction. 
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